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ABSTRACT 


A  theoretical  model  is  developed  based  on  an  iterated  perfect  gas 
inviscid-viscous  flow  field  which  includes  first-order  displacement 
{viscous  interaction),  transverse  curvature,  wall  slip,  and  temperature 
jump  in  addition  to  mass  transfer  effects.  The  effects  of  inviscid 
(tangent  cone)  and  viscous  (nonsimilar  laminar  boundary  layer)  flow 
field  matching  conditions  are  also  considered.  Numerical  results  are 
compared  with  experimental  and  analytical  results  of  King  and  Talbot 
(AIAA  J.,  1964)  for  a  5 -deg  half-angle  cone  at  M,,  =  3.  93  and  5.  64. 

The  predicted  viscous -induced  pressure  without  mass  transfer  and  the 
zero-lift  drag  with  and  without  injection  were  in  agreement  with  the 
experimental  data  within  experimental  uncertainty.  The  theoretical 
model  was  also  used  to  predict  zero-lift  drag  of  a  9-deg  half-angle 
cone  at  »  9.  37  and  10  and  Re^/in.  -  400  to  45,  000  for  a  range  of  wall- 
to-stagnation  temperature  ratios.  Again,  in  general,  without  mass 
transfer  the  predictions  were  within  experimental  uncertainty.  The  in¬ 
ability  of  the  theoretical  model  to  adequately  treat  nonuniform  mass 
transfer  distributions  is  discussed.  At  the  lowest  Reynolds  number  the 
effects  of  slip  were  most  significant.  At  all  conditions  the  effects  of 
inviscid-viscous  flow  field  matching  were  significant.  Experimental 
zero  injection  equilibrium  wall  temperature  distributions  and  cool-wall 
pressure  data  are  given  at  ~  10  and  Re^/in.  =  400  to  2600, 
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SECTION  I 
INTRODUCTION 


For  several  years  viscous  effects  on  sharp  and  blunt  slender  cones 
have  been  studied  in  the  Arnold  Engineering  Development  Center, 
von  Karman  Gas  Dynamics  Facility  (Refs.  1-3).  The  large  viscous- 
induced  drag  increment  at  zero-lift  observed  in  these  studies  encouraged 
Whitfield  and  Griffith  (Ref.  2)  to  consider  the  effects  of  mass  transfer  to 
reduce  the  drag.  They  experimentally  investigated  the  effects  of  tangen¬ 
tial  injection  of  gases  on  zero-lift  drag  at  =  10  and  21.  The  experi¬ 
mental  investigations  have  recently  been  extended  by  Little  and  Griffith 
(Ref.  4)  to  include  the  effects  of  normal  injection  on  zero-lift  drag  at 
=  10.  The  present  report  presents  an  analysis  of  some  of  the  latter 
experimental  results  after  developing  a  theoretical  model  and  making 
comparisons  with  some  experimental  and  theoretical  data  at  supersonic 
conditions. 

King  and  Talbot  (Refs.  5  and  6)  studied  the  effects  of  normal  gas  in¬ 
jection  on  the  pressure  distributions  and  zero-lift  drag  of  a  sharp  5-deg 
half-angle  cone  in  the  Berkeley  Low  Density  Wind  Tunnel.  Experimental 
data  were  taken  at  near  adiabatic  wall  conditions  (Tw/To  =  0.96)  and 
Re^/in.  •=  1765  and  6286  at  =  3.93  and  5,64,  respectively.  They  gave 
an  analysis  based  on  approximate  methods  which  is  described  in  more 
detail  below.  The  experimental  and  theoretical  data  of  King  and  Talbot 
are  compared  with  the  results  from  methods  developed  in  the  present 
report. 

King  and  Talbot  (Ref.  6)  and  Little  and  Griffith  (Ref.  4)  obtained 
experimental  data  for  binary  boundary  layers  with  injected  gases  other 
than  the  test  gases  (air  and  nitrogen)  used  in  the  wind  tunnel  experiments. 
For  example,  the  latter  investigators  also  used  helium,  argon,  and 
carbon  dioxide  as  injected  gases.  In  the  present  report,  however,  only 
the  cases  of  injection  of  air  into  air  or  nitrogen  into  air  or  nitrogen  will 
be  considered. 

The  experimental  drag  data  of  Little  and  Griffith  were  taken  with  a 
partially  porous  9-deg  half-angle  cone  at  =  9.  37  and  10  and  for  a 
span  of  Reynolds  numbers  Re^  /in.  =  388  to  42,  640.  The  estimated  range 
of  wall-to-stagnation  temperature  ratios  was  from  0.  2  to  0.  75.  Quite 
recently  some  experimental  data  were  obtained  on  the  equilibrium  wall 
temperature  variation  over  a  hollow  cone  in  the  AEDC-VKF  low  density 
tunnel  (Gas  Dynamic  Wind  Tunnel,  Hypersonic  (L)),  and  the  pressure 
was  measured  on  a  sharp,  water-cooled  cone  at  Re^/in.  =  388,  1602, 
and  2606  in  the  same  wind  tunnel.  These  recent  experimental  data  are 
also  included  in  this  report. 
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Under  certain  experimental  conditions  the  effects  of  slip  and  tem¬ 
perature  jump  were  significant,  and  the  present  analysis  presents  an 
approximate  treatment  of  those  effects  on  pressure  and  heat-transfer 
distributions  and  zero-lift  drag.  The  effects  of  slip  were  considered 
significant  upon  comparison  of  the  experimental  results  with  the  numer¬ 
ical  results  of  the  present  investigation  when  slip  and  temperature  jump 
effects  were  neglected. 

Before  further  considering  the  theoretical  model  used  in  the  present 
report,  the  experimental  conditions  which  were  studied  will  be  presented. 


SECTION  II 

EXPERIMENTAL  CONDITIONS  STUDIED 


2.1  WIND  TUNNELS 

The  experimental  data  used  in  the  present  study  were  obtained  at 
Berkeley  by  King  and  Talbot  (Refs.  5  and  6)  and  at  the  AEDC-VKF  by 
Little  and  Griffith  (Ref,  4).  Some  details  concerning  their  experimental 
investigations  are  given  in  this  section  for  the  convenience  of  the  reader. 

Three  wind  tunnels  were  used  in  the  experimental  studies:  viz.  the 
low  density  wind  ttinnel  at  Berkeley  and  two  AEDC-VKF  Mach  10  tunnels 
(Gas  Dynamic  Wind  Tunnels,  Hypersonic  (C)  and  (L)).  The  Berkeley 
tunnel  is  an  unheated  air,  continuous  flow,  free  jet  facility  using  axisym- 
metric  contoured  nozzles  to  produce  Mach  numbers  of  3.  93  and  5,  64. 
Tunnel  C  is  a  continuous  flow,  hypersonic  wind  tunnel.  Air,  heated  to  a 
temperature  sufficient  to  avoid  liquefaction,  expands  through  a  con¬ 
toured,  axisymmetric  nozzle  to  the  50-in,  -diameter  test  section.  The 
free-stream  Reynolds  number  was  varied  by  changing  the  free-stream 
stagnation  (reservoir)  pressure.  Tunnel  L  is  an  arc-heated,  continuous 
flow  facility  using  interchangeable  axisymmetric  contoured  nozzles  to 
produce  nominal  Mach  10  flow.  Nitrogen  was  used  as  the  test  gas. 

The  flow  conditions  from  the  experimental  investigations  which  are 
considered  in  the  present  report  are  given  in  Table  1.  The  experi¬ 
mental  conditions  are  given  in  the  first  row,  and  the  data  used  in  the 
niimerical  calculations  are  given  in  the  second  row.  For  most  cases 
the  perfect  gas  data  for  the  numerical  calculations  were  based  on  the 
measured  free-stream  pitot  pressure  and  (equilibrium)  free-stream 
stagnation  temperature.  The  equilibrium  stagnation  enthalpy  was  then 
determined,  and  the  ideal  gas  (y  =  1.4)  stagnation  temperature  was  com¬ 
puted.  Ideal  gas,  normal  shock,  and  one -dimensional,  isentropic  flow 
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relations  were  then  used  to  obtain  the  free -stream  and  reservoir  con¬ 
ditions.  This  procedure  gave  an  ideal  gas  stagnation  temperature 
somewhat  higher  than  the  corresponding  equilibrium  value,  and  this 
also  produced  other  differences  as  shown  in  Table  I. 


2.2  MODELS 

The  force  models  used  for  the  experimental  studies  were  sharp 
cones  with  partially  porous  surfaces.  King's  models  were  5 -deg  half- 
angle  cones  with  small  solid  base  and  tip  regions  (Fig.  la).  The  porous 
material  was  sintered,  powdered.  Monel*  metal  particles.  An  experi¬ 
mentally  determined  mass  flux  distribution  is  shown  in  Fig.  2.  King 
states  that  the  distribution  of  mass  flux  represents  a  circumferential 
average  and  should  be  regarded  as  qualitative  rather  than  quantitative. 
Surface  pressure  data  were  taken  with  mass  injection,  but  since  the 
porous  surface  near  the  taps  was  distorted  the  pressure  data  were 
affected.  Therefore,  the  only  experimental  pressure  data  of  King  and 
Talbot  used  in  the  present  report  were  taken  on  solid  (no  injection)  cone 
models . 

The  AEDC-VKF  force  model  was  a  sharp,  9 -deg  half-angle  cone 
with  a  base  diameter  of  0.  75 -in.  (Fig.  lb).  The  solid  nose  was  followed 
by  a  porous  afterbody  machined  from  General  Electric  50-percent  dense 
nickel  Foametal®.  The  impervious  nose  extended  over  0.57  of  the  (slant) 
length  of  the  cone.  Over  the  remainder  of  the  cone,  however,  the  con¬ 
stant  density  and  thickness  of  the  wall  gave  a  mass  flux  distribution 
constant  within  ±15  percent  as  obtained  from  a  hot-wire  anemometer 
survey.  Figure  2  also  shows  a  typical  mass  flux  distribution  for  the 
AEDC-VKF  model. 

The  zero  injection  drag  data  were  taken  after  the  model  had  been 
in  the  flow  a  sufficient  time  to  reach  approximately  90  percent  of  the 
equilibrium  wall  temperature  (as  indicated  from  changes  in  measured 
drag  data).  Although  no  experimental  wall  temperature  data  exist  for 
the  force  model,  the  values  quoted  below  for  a  hollow  stainless  steel 
model  are  believed  representative  for  the  zero-injection  drag  data  as 
well. 

The  pressure  was  measured  at  two  locations  on  a  water-cooled 
copper  model  (Pig.  Ic)  in  Tunnel  L.  The  thin  copper  wall  was  cooled 
to  the  water  temperature  {-  293“K),  and  the  temperature  was  assumed  to 
remain  constant  over  the  entire  surface. 

The  wall  temperature  was  measured  with  a  thin -wall  (t  <  0.010  in.) 
stainless  steel  model  as  shown  in  Fig.  Id.  The  output  of  thermocouples 
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welded  to  the  inside  surface  of  the  model  was  recorded  as  a  function  of 
time.  Wall  temperatures  reached  90  percent  of  the  equilibrium  value 
within  one  to  three  minutes,  and  equilibrium  was  reached  within  about 
ten  minutes. 


2.3  EXPERIMENTAL  DATA  AND  PRECISION 

The  primary  test  measurements  were  axial  force  and  surface  pres¬ 
sures  in  the  Berkeley  experiments  and  axial  force  with  and  without  mass 
transfer,  and  surface  pressures  and  wall  temperatures  without  mass 
transfer  in  the  AEDC-VKF  experiments.  The  estimated  precision  of 
these  data  is  shown  in  the  table  below: 


Source 

Force,  percent 

Pressure,  percent 

Pressure,  percent 

ITI  -  0 

m  >  0 

Berkeley 

±5 

±2 

±5 

Tunnel  C 

±7 

— 

— 

Tunnel  L 

±7 

•  » 

The  experimental  value  of  Co  for  each  value  of  injected  mass  flow 
used  in  the  numerical  calculations  is  shown  in  Table  II,  the  experimental 
pressure  data  are  shown  in  Table  III,  and  wall  temperature  data  are 
shown  in  Table  IV.  It  should  be  noted  that  the  drag  data  presented  do  not 
necessarily  correspond  to  actual  data  points  but  rather  were  taken  from 
fairings  of  the  data. 

Wall- to -stagnation  temperatures  are  also  tabulated.  For  the 
Berkeley  tests  the  data  were  from  readings  of  thermocouples  located  at 
the  base  of  the  model.  The  ratios  for  the  AEDC-VKF  zero  injection  data 
in  Tunnel  L  are  shown  in  Fig.  3,  and  the  data  from  Tunnel  C  are  from 
the  previous  data  of  Whitfield  and  Griffith.  The  accuracy  of  the  thermo¬ 
couple  output  is  about  ±2  percent;  however  no  attempt  has  been  made  to 
assess  effects  of  conduction  losses  through  the  wall.  With  injection  the 
assumption  is  made  that  the  cooled  wall  was  maintained  at  the  injection 
gas  temperature  (300“K). 

No  precision  is  quoted  here  for  the  pressure  data  in  Tunnel  L.  How¬ 
ever,  the  repeatability  of  the  measured  pressures  was  ±2  percent.  (Also, 
the  free-stream  pitot  pressure  for  the  conditions  -  9.37  and  Re^/in.  = 
1602  was  repeated  and  found  to  agree  with  previous  tunnel  calibrations 
within  1  percent. )  Relatively  large  corrections  to  the  measured  pres¬ 
sure  data  are  required  to  account  for  orifice  size  and  heat -transfer  effects. 
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and  these  effects  will  be  discussed  below.  Because  of  the  necessary 
corrections,  a  precision  for  the  pressure  data  was  more  difficult  to 
determine. 


SECTION  III 

SOME  THEORETICAL  CONSIDERATIONS 


3.1  KING  AND  TALBOT'S  ANALYSIS 

The  Berkeley  experimental  data  were  analyzed  by  King  and  Talbot 
using  (1)  the  similar  boundary -layer  solutions  of  Low  (Ref.  7),  (2)  the 
effects  of  transverse  curvature  reported  by  Yasuhara  (Ref.  8),  (3)  the 
inviscid  flow  field  from  the  approximate  tangent  cone  theory,  and 
(4)  the  matching  of  the  inviscid -viscous  flow  fields  based  on  d  S*/dx. 

The  "second-order"  boundary-layer  effects  were  assumed  to  be  un¬ 
coupled,  and  a  linear  combination  of  these  effects  was  added  to  the 
effects  of  mass  transfer  from  Low's  solutions.*  Also,  the  effects  of 
the  impermeable  nose  were  neglected  in  their  calculations. 

At  the  time  the  analysis  was  made  by  King  and  Talbot  (1962),  the 
above-mentioned  theories  and  procedures  were  perhaps  the  best  in 
existence.  However,  since  that  time  certain  theoretical  developments 
have  occurred  which  make  a  re-examination  of  their  experimental  data 
and  comparisons  with  their  analytical  results  of  interest.  This  seems 
especially  significant  since  their  experimental  data  were  taken  at 
supersonic  Mach  numbers  and  moderate  Reynolds  number  where  other 
second-order  boimdary-layer  effects  (viz.,  external  vorticity,  slip, 
and  temperature  jump)  should  be  small.  Therefore,  both  their  experi¬ 
mental  data  and  analysis  provide  interesting  data  to  test  new  theoretical 
models  designed  for  application  at  hypersonic,  cold-wall,  low  Reynolds 
number  conditions. 


3.2  PRESENT  ANALYSIS 

The  basic  theoretical  model  used  in  the  numerical  calculations  was 
an  iterated  perfect  gas  {y  =  1.4)  inviscid-viscous  flow  field.  The  viscous 
flow  field  was  computed  with  the  nonsimilar,  laminar  boundary -layer 
theory  of  Clutter  and  Smith  (Ref.  9)  which  includes  approximate  trans¬ 
verse  curvature  terms.  In  addition,  all  boundary-layer  calculations 


*Low's  flat-plate  similarity  solutions  were  transformed  with  the 
Mangier  factor  (3)^^^. 
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were  based  on  Sutherland's  viscosity  law  and  Prandtl  number  (Pr)  =  0.71= 
constant.  The  inviscid  flow  field  was  computed  from  the  approximate 
tangent -cone  theory.  The  matching  of  the  inviscid  and  viscous  flow 
fields  was  treated  in  two  ways:  (1)  using  the  usual  matching  condition 
Sell  -  dc  +  tan"‘  (dS*/dx)  and  (2)  (9eff  =  0c  +  tan~‘  (Ve/ug)  ,  where  Ve/uo 
was  computed  from  the  formula  proposed  by  Li  and  Gross  (Ref.  10).* 
Some  comments  on  the  matching  conditions  are  given  in  Appendix  I. 

The  inviscid -viscous  flow  field  was  iterated  until  there  was  negli¬ 
gible  chsinge  in  5*  (x)  or  (x) .  At  high  Reynolds  numbers  only  two 
iterations  or  three  boundary -layer  calculations  were  required.  At  low 
Reynolds  numbers  and  especially  at  high  injection  rates,  five  iterations 
or  six  boundary -layer  calculations  were  required.  Moreover,  signifi¬ 
cant  effect  of  matching  conditions  was  noted  on  the  number  of  iterations 
required  (the  number  being  less  with  ve/ug  since  with  this  matching  con¬ 
dition  the  inviscid  outer  flow  was  displaced  less  from  the  geometric 
wall  6c  (x)  =  constant).  The  method  used  for  iterating  the  flow  fields  is 
discussed  in  Appendix  II. 

It  should  be  noted  that  King  and  Talbot  used  a  graphical  method  for 
matching  the  inviscid  outer  flow  (tangent  cone)  with  the  viscous, 
boundary -layer  edge  conditions  (Low's  similar  solutions).  The  match¬ 
ing  was  done  simultaneously  step  by  step  along  the  body  rather  than 
obtaining  successive  approximations  to  the  entire  "effective"  body  as 
was  done  in  the  present  investigation, 

Because  of  the  difficulties  associated  with  iterating  the  inviscid- 
viscous  flow  field,  most  calculations  presented  in  this  report  were  based 
on  a  uniform  injection  (constant  blowing)  rate  rather  than  constant 
fw  -  lA'^*  as  required  for  similarity  (see,  e.g. ,  Low,  Ref.  7).  Uniform 
injection  more  nearly  represented  the  actual  mass  flux  distributions 
over  the  porous  region  of  the  model  in  the  =  10  experiments  and 
appears  to  be  a  reasonable  approximation  for  the  experimental  distribu¬ 
tions  of  King  and  Talbot.  The  effects  of  the  impermeable  nose  region 
were  investigated  for  the  King  and  Talbot  conditions  and  for  one  set  of 
conditions  at  =  10. 

The  experimental  data  of  King  and  Talbot  indicated  a  near-constant, 
wall -to -stagnation  temperature  ratio  Tw /To  =  0.96.  Most  of  the  numerical 
calculations  for  the  present  report  at  -  10  were  completed  before  the 
experimental  wall  temperature  data  were  available.  The  wall-to- 
stagnation  temperature  ratios  assumed  in  the  calculations  are  indicated 


♦Hereafter  the  matching  conditions  will  be  denoted  as  d  fiVA  or  vg/ue  . 
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in  the  figures  where  comparisons  with  experimental  data  are  given.  In 
all  calculations  the  wall  temperature  was  assumed  to  be  constant,  and 
the  effects  of  wall  temperature  distributions  such  as  shown  in  Fig.  3 
have  not  been  investigated. 

The  effects  of  slip  and  temperature  jump  were  investigated  at  all 
conditions.  The  boundary-layer  theory  of  Clutter  and  Smith  (Ref,  9) 
was  modified  to  include  these  effects.  For  more  direct  comparison 
with  others,  the  mean  free  path  was  evaluated  at  the  wall  temperature, 
Tvi,  rather  than  at  the  slip  condition  T(ij  =0)  ;  however,  this  is  not  a 
necessary  condition  for  the  numerical  procedure.  Details  of  the  pro¬ 
cedure  used  in  the  present  investigation  are  given  in  Appendix  III. 

It  should  be  noted  that  the  so-called  second-order  boundary-layer 
effects  have  all  been  treated  as  first-order  effects.  The  first-order 
momentum  and  energy  equations  included  (approximate)*  transverse 
curvature  terms,  and  the  boundary  conditions  were  modified  to  include 
slip  and  temperature  jump.  Moreover,  the  displacement  effect  is  con¬ 
sidered  a  first-order  effect  since  the  inviscid  and  viscous  flow  fields 
were  iterated  to  "convergence".  This  is  not,  however,  internally  con¬ 
sistent  since  the  outer  flow  was  treated  as  inviscid  (see  Van  Dyke, 

Ref.  1),  and  low  density  effects  on  the  shock  wave  (e.g. ,  finite  thick¬ 
ness,  curvature,  and  slip)  were  neglected  (see  Pan  and  Probstein, 

Ref.  12).  These  deficiencies  can,  in  principle,  be  eliminated  by  a 
proper  treatment  of  the  "outer"  flow. 


SECTION  IV  ' 
RESULTS  AND  DISCUSSION 


The  theoretical  model  used  in  the  present  report  was  tested  against 
King  and  Talbot's  experimental  data  and  analytical  results  at  =  3.  93 
and  5.  64.  The  theoretical  model  was  then  used  to  analyze  the  AEDC-VKF 
experimental  data  of  Little  and  Griffith  at  =  9,  37  and  10,  A  large 
volume  of  interesting  results  were  generated;  however,  the  data  pre¬ 
sented  herein  will  be  confined  to  some  of  the  more  significant  results. 


4.1  BERKELEY  CONDITIONS 

The  inviscid-viscous'  flow  field  was  iterated  as  described  above  and 
in  Appendix  II.  For  convenience  of  comparison  with  the  analyticad  results 


*See  the  revised  version  of  Clutter  and  Smith  (Ref.  9), 
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of  King  and  Talbot  at  =  5.  64  and  Re^/in.  =  6200,  the  ratio  6*/t„  was 
computed  without  mass  transfer  and  is  shown  in  Fig.  4.  Excellent 
agreement  was  found  with  King  and  Talbot's  results  using  the  d  5*/dx 
matching  condition  and  the  similar  solutions  of  Low.  The  matching  con¬ 
dition  ve/ue  predicted  a  slightly  larger  5*  since  the  inviscid  outer  flow 
was  displaced  less  by  the  boundary  layer,  and  thus  the  viscous  -induced 
pressure  was  less  which  in  turn  produced  larger  boundary -layer  and 
displacement  thicknesses.  Under  these  conditions,  transverse  curva¬ 
ture  had  a  negligible  effect  on  S*  (x) . 

The  viscous -induced  pressure  increment  (p  -  Pi„vVPinv  =  3.93 

and  5.  64  is  shown  in  Fig.  5.  The  effects  of  matching  condition  and 
transverse  curvature  <TVC)  were  also  considered.  Using  the  matching 
condition  ve/u*  gave  results  which  were  in  better  agreement  with  the 
experimental  pressure  data  of  King  and  Talbot,  especially  at  =  5.  64. 
The  effects  of  slip  and  temperature  jump  are  also  shown.  As  suggested 
by  King  (Ref.  5),  the  effects  of  slip  and  temperature  jump  on  the  pres¬ 
sure  distribution  were  indeed  small.  However,  we  will  see  later  that 
the  effects  on  the  zero-lift  drag  are  somewhat  larger. 

The  effects  of  uniform  mass  transfer  on  the  viscous-induced  pres¬ 
sure  increment  at  -  5.  64  are  shown  in  Fig.  6.  At  the  highest  injec¬ 
tion  rate  considered  (0.  09375  Ibm/hr),  the  effects  of  matching  condi¬ 
tions  are  also  shown.  At  high  injection  rates  the  effects  of  matching 
conditions  on  the  viscous -induced  pressure  increment  became  more 
important  especially  over  the  forward  portion  of  the  cone  as  shown  in 
Fig.  6. 

The  ratio  of  local  skin-friction  coefficients  Cf/Cf**'  at  =  5.  64  is 
shown  in  Fig.  7,  where  is  the  value  without  viscous  interaction, 
transverse  curvature,  or  mass  transfer.  The  separate  and  combined 
effects  of  viscous  interaction  and  transverse  curvature  at  zero  injection 
are  in  contrast  with  the  King  and  Talbot  results.  As  noted  by  King 
(Ref.  5),  should  be  greater  than  unity;  however, his  analysis 

predicted  <  1,  a  result  which  he  attributed  to  his  neglecting 

the  viscous -induced  axial  pressure  gradient  by  using  the  flat -plate  re¬ 
sults  of  Low  (Ref.  7).  The  effect  of  transverse  curvature  on  Cf  based 
on  the  approximate  theory  of  Clutter  and  Smith  was  in  substantial  agree¬ 
ment  with  King  and  Talbot's  prediction  using  the  theory  of  Yasuhara. 

The  coupled  and  uncoupled  effects  of  viscous  interaction  and  transverse 
curvature  onCf  could  be  seen  by  comparing  the  results  denoted 
VI  +  TVC^*’  with  (VI  +  TVC)'*^  The  comparison  (not  shown  in  Fig.  7) 
indicated  excellent  agreement.  This  tends  to  support  the  usual  assump¬ 
tion  of  linear  independence  of  these  second-order  effects.  The  results 
denoted  (VI  +  TVC)**\  however,  have  treated  both  as  coupled  first-order 
effects. 
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Three  injection  laws,  viz.,  <1)  ■  constant,  (2)  f„  -  constant, 

and  (3)  pwvv  -  from  the  experimentally  measured  distribution  of 
King  and  Talbot  at  =  5.64,  and  their  effects  on  friction  drag  are 
shown  in  Fig.  8.  The  effects  of  viscous  interaction  and  transverse 
curvature  coupled  with  the  effects  of  injection  laws  on  are  also 
shown.  One  significant  result  for  the  conditions  considered  was  that 
the  effects  of  the  impervious  nose  and  various  injection  laws  on  the 
friction  drag  were  negligible.  It  is  also  significant  that  the  prediction 
including  viscous  interaction  and  transverse  curvature  was  in  good 
agreement  with  the  faired  experimental  data  of  King  and  Talbot.  The 
latter  data  were  obtained  by  subtracting  the  integrated  experimental 
pressure  drag  from  the  experimentally  measured  total  drag  with  a  cor¬ 
rection  for  mass  injection  (cf. ,  King,  Ref,  5), 

Figure  9  shows  the  drag  components  and  total  drag  coefficient 
Cd  =  Cof  +  Cop  from  the  present  investigation  compared  with  the  experi¬ 
mental  data  and  analytical  results  of  King  and  Talbot  at  =  3.  93  and 
5.64.  •  The  effects  of  viscous  interaction,  transverse  curvature,  and 
matching  conditions  with  mass  transfer  are  shown  in  the  figure.  The 
present  results  are  in  substantially  better  agreement  with  the  experi¬ 
mental  data  than  were  the  predictions  of  King  and  Talbot.  In  fact,  the 
present  results  iterated  to  convergence  and  based  on  ve/ue  matching 
conditions  and  with  or  without  injection  predicted  total  drag  in  agree¬ 
ment  with  the  experimental  data  at  =  5.  64  within  the  experimental 
uncertainty.  This  encouraging  result  should  be  tempered  with  the  ob¬ 
servation  that  the  present  analysis  cannot  be  extended  to  much  higher 
injection  rates.  Without  viscous  interaction,  "blow-off”  occurred  far 
forward  on  the  body  at  m  -  0.12  Ibm/hr,  whereas  King's  experimental 
data  extended  to  m  =  0,  32  Ibm/hr  before  injection  caused  the  total  drag 
to  reach  a  minimum. 

The  effects  of  slip  and  temperature  jump  on  the  zero  injection  drag 
at  =  3.  93  and  5.  64  are  shown  in  Fig.  9  (the  procedure  used  is  dis¬ 
cussed  in  more  detail  below  for  =  10  conditions).  At  both  Mach  num¬ 
bers,  the  effects  of  slip  and  temperature  jump  were  about  5  percent,  and 
this  should  be  contrasted  with  the  almost  negligible  effect  on  the  viscous - 
induced  pressure  increment  (see  Fig.  5).  If  the  present  theoretical 
model  is  appropriate  with  respect  to  slip  and  temperature  jump,  it 
appears  once  again  that  experimental  and  theoretical  zero-lift  drag 
studies  might  be  a  useful  tool  for  investigating  another  second-order 
boundary-layer  effect.  Of  course,  much  work  remains  to  solve  the 
leading  edge  problem  as  well  as  the  kinetic  theory  problems,  such  as 
reflection  and  accommodation  coefficients.  Although  the  predicted  effects 
of  slip  at  the  Berkeley  conditions  were  approximately  equal  to  experi¬ 
mental  uncertainty,  the  effects  may  be  large  enough  to  permit  experi¬ 
mental  verification  at  similar  conditions. 
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4.2  AEDC’YKF  CONDITIONS 

The  theoretical  model  was  applied  to  the  experimental  data  of  Little 
and  Griffith  at  =  10.  The  effects  of  iteration  on  the  displacement 
thickness  and  corresponding  viscous -induced  pressure  increment  with 
mass  transfer  at  Re^/in.  =  396  are  shown  in  Fig.  10.  With  injection  the 
displacement  thickness  was  substantially  increased,  and  this  in  turn 
produced  a  large  variation  in  the  viscous -induced  pressure  increment 
with  iteration.  It  should  be  noted  that  the  solution  had  not  converged 
with  five  iterations;  however,  the  converged  solution  was  bounded  by  the 
fourth  and  fifth  iterations. 

Figure  11  shows  the  effects  of  different  injection  laws  on  the  zeroth 
iteration  values  of  S  and  S*.  The  difficulty  of  iterating  the  inviscid- 
viscous  flow  field  near  the  beginning  of  the  injection  region  (x  =  0.  57  to 
0.  75)  prevented  us  from  considering  further  the  experimental  distribu¬ 
tion.  Substantial  differences  in  d  and  S*  caused  by  the  two  distributions 
(experimental  and  constant  mass  transfer)  are  evident  in  the  figure. 

Figure  12  shows  the  effects  of  mass  transfer  on  the  displacement 
and  total  boundary -layer  thicknesses  at  Re^/in.  =  396.  The  thicknesses 
were,  of  course,  increased;  however,  the  character  of  the  thicknesses 
was  not  substantially  changed  at  this  moderately  high  injection  rate. 


4.3  EFFECTS  OF  SLIP  AND  TEMPERATURE  JUMP 

Before  presenting  the  comparisons  of  the  numerical  results  with  the 
experimental  data  for  the  effects  of  mass  transfer  on  the  drag  coefficient 
at  =  9.37  and  10,  the  effects  of  slip  and  temperature  jump  on  the 
numerical  results  are  considered. 

Figure  13  shows  the  ratio  of  the  wall  slip  velocity  to  free-stream 
velocity  and  temperature  jump  at  =  10  and  the  lowest  Reynolds  num¬ 
ber  for  which  experimental  data  were  available.  The  effects  of  match¬ 
ing  conditions  are  shown,  and  the  effects  of  viscous  interaction  are 
implicit  since  the  calculations  with  slip  and  temperature  jump  were 
started  from  the  converged  no-slip  solutions  (see  Appendix  III).  The 
temperature  jump  at  the  wall  especially  deserves  a  few  comments. 

Near  x  =  o  where  the  slip  velocity  is  largest,  the  temperature  jump  is 
smallest  -  a  surprising  result.  However,  the  temperature  jump  is 
dominated  by  the  quantity  (db/d7/)^^^  (see  Appendix  III).  From  the  energy 
equation,  dh/ dy  =  Heg'  -  Ue*  F  t"  ,  where  prime  denotes  differentiation 
with  respect  to  rf .  Thus  with  slip  there  exist  values  of  f',  g'  and  f'"such 
that  dh/dtj  =  0  or  h  ()j  =  0)  =  Kk  .  Also,  the  solution  based  on  dSVdx 
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matching  conditions  reached  a  maximum  temperature  jump  near  x  =  0.7 
and  began  to  decrease.  Of  course,  the  asymptotic  limit  T(i?  =  0)  =  Tw 
must  be  reached  as  x  « j  however,  for  this- low  Reynolds  number  con¬ 
dition  there  were  large  slip  and  temperature  jump  effects  over  the 
entire  cone. 

The  effects  of  slip*  and  temperature  jump  on  S*  and  S  are  shown  in 
Fig.  14.  Also,  the  effects  of  mass  transfer  and  slip  on  the  pressure 
distribution  are  shown  in  Fig.  15.  At  zero  injection,  the  effects  of  slip 
on  the  pressure  distribution  became  negligibly  small  as  x  ^  1 , 

The  skin-friction  coefficient,  Cf^,  and  Stanton  niunber,  St^,  are 
shown  in  Figs.  16a  and  b,  respectively.  The  results  with  mass  transfer 
and  slip  are  shown,  and  the  effects  of  viscous  interaction  and  transverse 
curvature  were  included  in  the  calculations.  In  accord  with  the  pre¬ 
viously  discussed  results  with  slip,  the  effects  of  Cr^  and  St^  were  sub¬ 
stantially  reduced  as  x  increased;  however,  it  is  interesting  to  note 
that  the  effects  of  slip  were  somewhat  greater  on  St^  than  on  Ci  ,  where- 
as  the  effects  of  mass  transfer  on  the  two  coefficients  were  nearly  equal. 

Comparisons  of  experimental  pressure  data  and  theoretical  pre¬ 
dictions  are  shown  in  Fig.  17,  The  uncorrected  experimental  data  are 
indicated  with  flagged  symbols.  The  data  were  corrected  for  orifice 
size  after  Potter,  Kinslow,  and  Boylan  (Ref.  13),  and  those  results  are 
also  shown.  The  predictions  from  the  present  investigations  are  shown, 
and  for  comparison  the  inviscid  sharp  cone  values  and  the  free -molecule 
values  after  Hayes  and  Probstein  (Ref.  14)  are  also  shown.  The  present 
iterated  flow  field  and  free-moleciile  flow  predictions  both  overestimate 
the  measured  surface  pressure.  As  noted  earlier,  the  orifice  correction 
is  large  and  the  experimental  uncertainty  is  currently  unknown.  The 
trends  of  both  the  uncorrected  and  corrected  surface  pressure  data  are 
as  expected,  i,  e. ,  the  pressure  tends  toward  the  free-molecule  value 
as  X  -  0.  The  agreement  between  experimental  data  and  numerical  re¬ 
sults  is  encouraging;  however,  further  work  is  required  to  better  define 
the  experimental  uncertainty  and  correction  factors  (orifice  size  and 
thermal  transpiration)  as  well  as  the  applicability  of  the  theoretical 
model  at  these  low  Reynolds  number  conditions. 

It  is  also  worth  noting  that  the  difference  between  the  inviscid  sharp 
cone  value  and  the  free-molecule  value  is  primarily  due  to  the  large 
wall-to-free-stream  temperature  ratio.  The  wall  temperature  effect 
increases  with  increasing  Tw/T„  (see  Hayes  and  Probstein,  Ref.  14). 


*Hereafter  slip  will  be  used  to  denote  both  slip  and  temperature 
jump  effects. 
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Comparisons  of  the  experimental  drag  data  at  ~  10  with  the  pre¬ 
dictions  from  the  present  investigation  are  shown  in  Fig,  18.  The 
effects  of  mass  transfer,  viscous  interaction,  transverse  curvature, 
matching  conditions,  slip,  and  temperature  jump  are  indicated.  Only 
the  final  converged  solutions  are  shown  except  at  the  lowest  Reynolds 
nximber.  The  results  of  each  iteration  and  the  combined  effects  of  mass 
tremsfer  and  slip  at  toe  lowest  Rejmolds  number  are  shown  in  Pig.  19  to 
indicate  the  large  changes  in  the  total  drag  between  iterations.  The 
effect  of  the  term  pw  vwr„/peUere  on  the  matching  condition  with  and  with¬ 
out  slip  is  also  shown  in  the  figtire.  The  effect  of  the  term  on  the  total 
drag  was  about  5  percent  and  in  a  direction  away  from  the  experimental 
results. 

It  is  interesting  to  note  that  for  the  same  total  mass  flow  rate,  m , 
the  effects  of  injection  on  the  drag  increase  with  decreasing  Reynolds 
number.  This  is  also  in  agreement  with  the  experimental  results  of 
Whitfield  and  Griffith  (Ref.  2)  for  the  case  of  tangential  injection.  It 
thus  appears  that  the  effects  of  mass  transfer  (either  normal  or  tangen¬ 
tial)  are  more  important  on  inviscid-viscous  flow  fields  under  low 
Reynolds  number  conditions. 

In  considering  the  comparison  of  niunerical  results  and  experi¬ 
mental  data  with  mass  transfer,  the  basic  deficiency  in  the  theoretical 
model  should  be  remembered.  Uniform  mass  transfer  was  assumed 
over  the  entire  cone,  whereas  in  the  AEDC-VKF  experiments  the  for¬ 
ward  32,5  percent  of  the  cone  surface  area  was  impervious  to  mass 
transfer.  Moreover  a  theoretical  model  with  impervious  nose  and 
porous  afterbody  would  probably  predict  an  increase  in  total  drag  be¬ 
cause  of  (1)  an  increase  in  pressure  and  friction  drag  over  the  im¬ 
pervious  nose  and  (2)  an  increase  in  pressure  drag  over  the  porous 
afterbody.  It  is  not  clear  whether  the  friction  drag  would  be  increased 
or  decreased  over  the  porous  region  since  the  increased  mass  flux 
would  tend  to  reduce  the  skin  friction;  however,  the  thicker  boundary 
layer  would  tend  to  increase  the  viscous  and  transverse -curvature - 
induced  friction  drag.  Therefore  considering  the  deficiency  in  the 
theoretical  model  with  regard  to  nonuniform  mass  transfer  distribu¬ 
tions,  the  results  of  the  present  predictions  are  considered  to  be  satis¬ 
factory. 

A  comparison  of  numerical  results  and  experimental  data  for  all 
the  zero-lift  drag  data  without  mass  transfer  is  shown  in  Fig.  20.  For 
comparison,  the  previously  published  experimental  data  of  Whitfield 
and  Griffith  are  also  shown.  The  drag  data  are  shown  against  the 
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V 

hypersonic  interaction  parameter  s  ’  proposed  by 

Whitfield  and  Griffith.*  In  general,  the  numerical  results  based  on  the 
matching  condition  ve/ue  and  including  slip  agree  with  the  experimental 
data  within  experimental  uncertainty. 


SECTION  V 

CONCLUDING  REMARKS 


When  the  theoretical  model  developed  in  the  present  report  was 
used  to  predict  the  viscous -induced  pressure  increment  and  zero-lift 
drag  of  a  sharp,  5 -deg  half-angle  cone  at  the  Berkeley  conditions 

=  3.  93  and  5.  64),  the  results  were  In  better  agreement  with  the 
experimental  data  than  were  the  previous  analytical  results  of  King  and 
Talbot  (Refs.  5  and  6).  This  was  especially  true  for  the  effects  of  mass 
transfer  on  the  total  drag  at  moderately  high  injection  rates. 

The  same  theoretical  model  was  used  to  predict  viscous  effects  on 
a  9-deg  half-angle  cone  at  =  9.  37  and  10  over  a  range  of  Reynolds 
numbers  and  wall-to-stagnation  temperature  ratios.  At  zero  injection, 
in  general  the  predictions  of  zero-lift  drag  agreed  with  the  experi¬ 
mental  data  within  experimental  uncertainty.  However  with  injection, 
the  basic  deficiency  of  the  theoretical  model,  viz,,  neglecting  the 
effects  of  a  large  impermeable  nose,  was  more  important.  The  effects 
of  inviscid -viscous  flow  field  matching  conditions,  slip,  and  tempera¬ 
ture  jump  as  applied  in  the  present  report  were  found  to  have  a  signifi¬ 
cant  effect  on  viscous -induced  pressure  and  friction  drag  increments. 
The  comparison  of  experimental  and  predicted  surface  pressures 
showed  substantial  disagreement,  and  the  need  for  further  investigation 
is  indicated.  Without  mass  transfer  the  theoretical  model  of  this  report 
is  applicable  to  arbitrary,  pointed,  two-dimensional,  and  axisymmetric 
bodies  as  long  as  boundary -layer  theory  is  applicable  and  an  inviscid 
outer  flow  exists  which  can  be  approximated  by  tangent  wedge  or  cone 
theory.  Under  the  rarefied  flow  conditions  at  =  10  considered  in 
this  report,  the  above -noted  requirements  might  not  have  been  met. 
However,  the  comparisons  shown  are  encouraging,  and  further  experi¬ 
ments  and  numerical  studies  are  required  to  better  define  the  effects  of 
wall-to-stagnation  temperature  ratio  and  wall  temperature  distributions 
and  the  nature  of  the  outer  flow  field  and  shock  wave. 


*Some  additional  comments  and  a  derivation  of  the  parameter 
are  given  by  Lewis  and  Whitfield  (Ref.  3). 
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APPENDIX  I 

INVISCID-VISCOUS  FLOW  FIELD  MATCHING  CONDITIONS 


The  effects  of  mass  transfer  on  inviscid-viscous  flow  field  matching 
conditions  have  been  considered  by  Mann  (Ref.  15)  and  Thyson  and 
Schurmaim  (Ref.  16).  More  recently,  Li  and  Gross  (Ref.  10)  considered 
the  combined  effects  of  mass  transfer  and  transverse  curvature  by  inte¬ 
grating  the  continuity  equation  from  y  =  0  to  y  =  S  (x) .  In  so  doing  they 
obtained  the  following  matching  conditions  for  a  sharp  cone  in  hyper¬ 
sonic  flow: 


Li  and  Gross  then  state  "...  In  hypersonic  boundary  layers,  S* 
then  Eq.  (I-l)  becomes 

_  Pw  *'y(  *1 6 

~  P#  "e  d* 

(which)  provides  the  necessary  condition  ..." 


(1-2) 


For  two  of  the  conditions  studied  in  the  present  report,  the  follow¬ 
ing  data  were  obtained  from  the  zeroth  iteration  without  mass  transfer 
at  St  =0.  90: 

^C.deg  Tw/To  S*/S  tBn^‘ (d8*/<lx),  deg  lan“*  (ve/ue).  deg 
5.64  6200  5  0.96  0.78  1.6B  1.01 

10.0  396  9  0.20  0.71  4.12  1.97 


The  effect  of  the  terms  neglected  by  Li  and  Gross  under  the  hypersonic 
approximation  caused  a  40-  and  50-percent  change  in  tan~*  (d8*/dx)  at 
=  5.64  and  10,  respectively.  This  corresponds  to  a  reduction  of  the 
viscous -induced  pressure  term  p/pinv  -  1  of  39  and  56  percent  at 
=  5.  64  and  10,  respectively.  It  therefore  appears  significant  and 
desirable  to  evaluate  the  complete  Eq.  (I-l)  rather  than  using  the  hyper¬ 
sonic  approximation  and  the  resulting  Eq.  (1-2). 

The  effect  of  the  term  Pw  v^rw/poUeCe  in  Eq.  (1-2)  was  investigated 
at  =  10  and  Re^/in.  =  396  where  the  effect  on  ve/ue  was  a  maximum. 
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Under  those  conditions,  the  effect  on  total  drag  was  about  5  percent. 
Therefore  for  more  direct  comparison  with  King  and  Talbot, 

Ve/ue  =  dSVdx  (1-3) 

was  used  in  most  of  the  calculations.  The  inviscid  flow  field  was  com¬ 
puted  from  the  approximate  tangent  cone  theory.  The  pressure  on  the 
effective  cone  of  half-angle  fleff  =  +  tan“‘(ve/ue)  was  determined 

according  to  either  Eq.  (I-l)  or  (1-3),  and  matching  is  denoted  ve/ue 
or  dsVdx  ,  respectively. 
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APPENDIX  II 
ITERATION  SCHEME 


In  order  to  consider  the  coupled  effects  of  the  various  second-order 
terms  including  viscous  interaction,  it  was  necessary  to  iterate  succes¬ 
sive  boundary-layer  solutions  until  a  converged  solution  was  obtained. 
Regardless  of  whether  the  matching  condition  was  based  on  dsVdx  or 
ve/uo  I  it  was  difficult  to  obtain  sufficiently  smooth  pressure  and  velocity 
distributions.  Three-point  numerical  differentiation  of  S  (x)  and  S*  (x) 
was  satisfactory  at  the  lower  Mach  numbers,  =  3.  93  and  5.  64,  but 
was  unsatisfactory  at  =  9.  37  and  10  since  small  variations  in  the 
smoothness  of  the  9eff  profile,  inherent  in  the  numerical  differentiation 
scheme,  caused  appreciable  oscillations  in  succeeding  boundary -layer 
solutions.  Several  schemes  for  fitting  B  (x)  and  S*  (x>  with  polynomials 
were  tried,  and  although  good  fits  of  the  profiles  were  obtained  the  cor¬ 
responding  analytical  derivatives  were  poorly  behaved. 

The  S  and  6*  calculated  data  were  fitted  with  a  power  law  and 
exponential  equation  of  the  form 

5  =  ax**  -  exp  (ex  -  d) 

This  equation  fitted  the  data  well  (average  error  about  ±0.  2  percent)  for 
the  majority  of  the  cases  considered.  Also,  the  resulting  derivative 

dS/dx  =  abx^“*  -  c  exp  (cx  -  d) 

was  well  behaved.  However,  this  form  of  equation  characteristically 
gave  a  poor  solution  near  x  =  0,  and  to  account  for  this,  6en  was 
graphically  extrapolated  over  the  region  0  ^  x  <  0.1 .  This  gave  reason¬ 
able  values  of  pressure  and  velocity  for  starting  the  solution  and  had,  a 
negligible  effect  on  the  total  zero-lift  drag  since  only  1  percent  of  the 
cone  surface  area  was  within  the  extrapolated  region. 

For  the  cases  with  mass  transfer  at  the  wall,  the  problem  was 
further  complicated.  Solutions  based  on  the  experimental  injection  laws 
(at  ~  10  injection  started  at  x  =  0.57)  proved  impossible  to  iterate  by 
the  method  of  successive  approximations  since  the  rapid  changes  in 
dS*/dx  encountered  near  the  start  of  the  injection  region  induced  serious 
oscillations  into  the  solution.  In  the  cases  which  were  iterated,  the 
mass  flux  was  computed  from  p^vy,  =  m/A„  ,  where  m  is  the  total  (experi¬ 
mental)  mass  transfer  rate  and  Aw  is  the  total  surface  area  of  the  cone. 
The  resulting  average  mass  flux  was  assumed  constant  over  the  entire 
surface  of  the  cone.  This  eliminated  the  oscillations  and  permitted 
inviscid -viscous  flow  field  iterations  to  be  made.  When  the  mass  flux 
reached  a  sufficiently  large  value,  boundary-layer  "blow-off"  (f  ^  -*  0) 
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occurred  on  the  cone.  It  was  possible,  however,  to  obtain  iterated  flow 
field  solutions  at  moderately  high  injection  rates  if  the  rapid  increase 
in  5  (x)  and  B*  (x)  as  x  -*  1  was  ignored.  This  is  justified  on  the  basis 
that  the  corresponding  viscous -induced  pressure  acts  to  reduce  a  rapid 
increase  in  S*  (x)  .  It  is  also  noted  that  each  calculation  is  to  be  con¬ 
sidered  as  an  iterate  or  approximate  solution  to  the  final  converged 
solution,  and  thus  an  exact  fit  of  5*  (x)  was  not  necessarily  the  most 
desirable  for  each  successive  approximation  at  moderately  high  blow¬ 
ing  rates. 

With  the  "effective"  body  obtained  as  described  above,  the  "effective" 
cone  angle  fleff  aJid  matching  condition  dS*/dx  or  ve/ue  were  used  with  the 
ideal  gas  iy  =  1.4)  tangent  cone  solutions  to  determine  pe  W  and  ue  (x) . 

The  boundary -layer  calculations  then  proceeded  in  the  usual  manner. 
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APPENDIX  III 

SLIP  AND  TEMPERATURE  JUMP  CONDITIONS 


The  effects  of  slip  and  temperature  jump  at  the  surface  of  blunt- 
nosed  bodies  have  been  studied  by  several  investigators  in  recent  years 
(see,  e.g..  Refs.  11  and  17).  Recently  Pan  and  Probstein  (Ref.  12) 
and  Li  and  Street  (Ref.  18)  have  considered  an  impermeable,  semi¬ 
infinite  Qat  plate  under  rarefied  flow  conditions.  The  results  of  the 
latter  investigations  varied  greatly  on  the  prediction  of  the  beginning  of 
the  region  of  continuum,  strong  interaction  plus  wall  slip  and  tempera¬ 
ture  jump.  Moreover,  most  authors  have  evaluated  the  mean  free  path 
at  the  wall  conditions,  say  Tw  and  pw,  rather  than  at  the  slip  conditions 
T  (t;  =  0)  and  Pw  .  This  usual  procedure  was  also  used  in  the  present 
report  but  is  not  a  necessary  condition  for  the  numerical  calculations. 


The  wall  slip  velocity  and  temperature  jump  were  computed  accord¬ 
ing  to  the  formulas  given  by  Street  (Ref.  19)  and  Van  Dyke, (Ref.  11): 

u(0)  =  (/xR/p)  [(;rT/2R)‘^  d^/dy  +  (3/4)  5T/dx]  (lll-l) 

T(0)  -  Tw  (fiR/p)  (15/8)  (;rT/2R)‘'^'  5T/<?y  (III-2) 

where  the  right-hand  sides  of  both  equations  were  evaluated  at  the  wall 
conditions,  the  hard  sphere  model  was  used  to  evaluate  constants  in 
both  formulas,  and  the  reflection  and  thermal  accommodation  coeffi¬ 
cients  were  assumed  to  be  unity. 


Under  the  transformations  used  by  Clutter  and  Smith  (Ref.  9), 

X  =  X  and  =  (ueL/p*u^x)^  p  dy  (III-3) 

Eqs.  (III-l)  and  (I1I-2)  become,  respectively, 

‘4 


(III- 4) 


=  0)  -  ^  -i[(i  ^  K  R...L  f  I  » 

+  ~  y  d  /  h  \  1 

4  y  -  1  u,  a*  (  '’ref/wj 

K,,  .  0)  -  h.  .  h.  R...  .  ^  ±f  ^ 


The  local  skin- friction  coefficient 


% 


Cf^  =  2  rw/p,u,’  =  2(ue/u,)’  (1/x)  (pwpw/p,/i^>  f  w 


21 


AEDC.TR-66-37 


is  not  altered  because  of  slip  and  temperature  jump;  however,  the  heat- 
transfer  rate  to  the  surface  becomes 

qw  =  [(ft/Pr)  d\i/dy  4  /xu  ((3u/f)y )  1  ^ 

and  with  the  Stanton  number  defined 

Stj  =  ilw/p,  (He  ~  hw) 


the  ratio  St,/Cf  becomes 


St, 

Cf 


^  iii-  F—  i'  (\  -  J_'ll 

(l-g»)  -J.  [Pr  r;  ^  He  Pr^J 


(III-6) 


The  Clutter-Smith  theory  (Ref.  9)  was  modified  to  include  Eqs.  (III-4), 
(III-5),  and  (III-6),  However,  before  sharp-nosed  body  calculations 
could  be  made  it  was  necessary  to  consider  the  effects  of  the  singularity 
at  X  =  0.  The  order  of  the  singularity  for  a  pointed  body  is  such  as  to 
preclude  internally  consistent  calculations  at  that  point.  However,  the 
boundary-layer  equations  are  not  applicable  at  x  =  0  and  a  small  region 
near  x  =  0  was  in  near  free-molecule  flow.  The  extent  of  the  noncon¬ 
tinuum  region  for  a  sharp  cone  in  rarefied  flow  has  not  been  investi¬ 
gated.  However,  Pan  and  Probstein  have  considered  the  problem  for  a 
sharp  flat  plate,  and  their  estimates  for  the  conditions  of  the  present 
investigation  are  of  interest. 


Pan  and  Probstein  gave  the  criterion  /M^  =  8  for  the  beginning 

of  the  "boundary  layer  plus  normal  pressure  gradient"  region.  This 
region  lies  downstream  of  the  near  free-molecule  and  transition  regions 
and  upstream  of  the  hypersonic  strong  interaction  region  (see  Ref.  12). 
The  region  was  treated  by  Pan  and  Probstein  as  continuum  flow  with 
wall  slip  and  temperature  jump,  and  they  investigated  the  effects  of  slip, 
curvature,  and  thickness  of  the  shock  wave  in  addition  to  the  effects  of 
the  normal  pressure  gradient.  (In  the  present  report  all  transport 
effects  on  the  shock  wave  and  the  normal  pressure  gradient  effects  were 
ignored. ) 


Using  the  Pan  and  Probstein  criterion, 

"•-■I  ' _ 

with  /M^  =  8  ,  for  the  Tunnel  L  conditions  the  initial  location  x-, 

was  computed  as  follows: 
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TUNNEL  L 


He^/in. 

T„/To 

in. 

xi 

9.37 

2869 

0.6 

0.0048 

7.0 

0.011 

9.37 

1632 

0.4 

0.0085 

5.S 

0.021 

10 

396 

0.2 

0.037 

6.8 

0.105 

One  thus  sees  that  except  at  the  lowest  Reynolds  number,  the  non¬ 
continuum  region  was  very  small  in  extent.  Pan  and  Probstein  also 
predicted  that  the  region  of  strong  viscous  interaction  plus  slip  should 
begin  at  Re„  ~  40  to  45.  Thus  at  the  lowest  Reynolds  number 

approximately  one-half  of  the  cone  was  within  the  (flat-plate)  noncon¬ 
tinuum  and  boundary -layer  plus  normal  pressure  gradient  regions.  For 
the  other  conditions  no  more  than  0. 125  of  the  length  or  0.  0156  of  the 
area  of  the  cone  was  upstream  of  the  (flat -plate)  strong  interaction 
region.  Also  without  slip  the  strong  interaction  region  of  a  cone  is 
much  less  than  for  a  flat  plate. 

Therefore  in  the  present  study  the  noncontinuum  and  boundary  layer 
plus  normal  pressure  gradient  regions  were  neglected.  The  boundary- 
layer  solutions  with  slip  and  temperature  jump  were  begun  as  near  x  =  0 
as  practical  (usually  at  x  -  0,  025)  using  the  modified  Clutter -Smith 
theory  and  the  tangent  cone  theory.  The  calculations  were  started  from 
the  converged  no-slip  solutions.  At  the  nonsimilar  terms  (e.  g. , 
dl'/dx,  dg/dx)  in  the  momentum  and  energy  equations  were  set  to  zero. 
However,  the  transverse  curvature  term,  which  in  the  original  Clutter - 
Smith  machine  program  was  also  set  to  zero  at  the  initial  station,  was 
retained  in  the  modified  theory.  The  boundary  layer  with  slip  and  tem.- 
perature  jump  was  also  iterated  as  the  no-slip  solutions  were,  and  the 
method  of  iteration  is  discussed  in  Appendix  11. 
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All  Dimensions  in  Inches 


««  Berkeley  Low  Density  Tunnel  Force  Model 


All  Dimensions  in  Inches 


b,  AEDC-VKF  Force  Model 
Fig.  1  Force,  Proisvre,  ond  Woll  Temperature  Models 


AEDC-TR-66<37 


Ftg.  1  Concluded 
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Fig.  2  Comparison  of  Gas  Injection  Laws 
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Experimental  Data 
AEDC-VKF  Tunnel  L 


Fig.  3  Equilibrium  Wali  T«mp»rature  Distributions  ov«r  a  Hollow  9-deg  Half-Angle  Cone 


AEDC-TR.66.37 


M(n-3.93 
Rejp/m.  « 17© 

T^y/TQ«0.96INoSlipl 
Xg,  ”0.0032  in. 
Zero  Injection 
Numerical  Results 


Sym' 

thatching 

Condition 

Iteration 

Source 

d6*/dx 

2 

Present  Results 

d6*/dx 

- 

King  and  Talbot' 

V“e 

2 

Present  Results 

V*Je 

2s  +  ls 

Present  Results 

TVC 


X 


a.  =  3.93 

Fig,  5  Effects  of  Viscous  Interaction  and  Transverse  Curvature 
on  Sharp  Cane  Pressure  Distribution 


30 


AEDC-TR-fi6-37 


■5.M 

Re^j/ln. 

OO 

■ 

t 

-0.038 

-0.96  (No  Slip! 

^OJ 

-0.0012 

Zero  Injection 

Numerical  Results 

Matching 


31 


T^o"0.96 
Sym  Source 


_  TVC  Iteration  -  Present  Results 


No 


No  Interaction 


I _ I _ I _ I _ I _ I _ 1 _ I _ I _ I _ I 

3  0.1  0.2  0,3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

X 


Fig.  7  Effects  of  Viscous  Interaction  and  Transverse  Curvature  on  Skin-Friction  Coefficient 


AEDC-TR-66-37 


Mgj  =  5. 64 
Re^/in.  *=  6200 

Numerical  Results 

Injection  Law 
Experimental  ‘  Constant 


PwVw "  Constant 


No  VI  orTVC 
VI  Only 
TVCOnly 
TVC  + VI 


0.24 


I 


0.2, 


Flags  Denote  Number  of  Iterations 
No  Slip 


'Df 


0.1 


a 


X 


-Faired  Experimental  Data 
after  King  and  Talbot 


o 

X 


a 

o 

+ 


J_ 


X 


0 

V 

Q 

Q 


0.02  0.04  0.06  0.08  0.10  0.12  0.14 

ifi,  Ibm/hr 


Fig.  8  Effects  of  Vorious  Injection  Lows  on  Friction  Drag  Coefficient 


34 


AEDC-TR-66-37 


M^-3.93 
Re^j/in.  ■  1765 

0.96  (No  Slip! 


Experimental  Data 

Berkeley  Low  Density  Tunnel 
King  and  Talbot 


Sym 


& 


Matching 

Condition 

d6*/dx 
<16 '/dx 

V“e 


Numerical  Results 


Injection 

Law 

•  Const. 

“  Const. 

Pv„Vy,  ■  Const. 

Flag  Denotes  Slip 


Source 

Present  Results 
King  and  Talbot 
Present  Results 


a.  -  3.93 
Fig.  9  Drqg  Components 


35 


AEDC-TR-66-37 


Re^j/in.  -6200 
T,/To-0.96<NoSlip) 


Experimental  Data 
Berkeley  Low  Density  Tunnel 


—  King  and  Talbot 

Numerical  Results 

Matching  Injection 

Condition  Law 

Source 

— 

dfi«/dx 

Pw'^w  ■  Const. 

Present  Results 

— 

d6*/dx 

■  Const. 

King  and  Talbot 

A 

Vi^e 

PvjV^,  •  Const. 

Present  Results 

10  dag 

■  J.  5  in.  ♦] 

Fiag  Denotes  Slip 

b.  =  5.64 
Fig.  9  Concluded 


36 


Re^j/in.  =396 

X(jj  =  0.037  In. 
rh  =  0. 144  Ibm/hr 
d6*/dx  Matching  Condition 
No  Slip 


Iteration 


1.50  0.60  0.70 


a.  On  Displacement  Thickn«ss 
Fig.  10  Eifects  of  Iteration  with  Mass  Transfer 


AEDC-TR-66-37 


Woo'/^-<>-326 

X3,=  0.037  In. 
m  =*0. 144  Ibm/hr 
(16 “/dx  Matching  Condition 
No  Slip 


Iteration 


II 


ZC*99*ai-Da3V 


AEDC.TR-66-37 


M^-10 
Reojj/ln.  -  396 

T^^,n-o"0.2(Nosnp) 
X^- 0.037  In. 
Zero  Injection 


Iteration 


0.326 


T^y^o-0.2{NoSlip) 


Fig.  14  Ef^cts  of  Slip  and  Temperature  Jump  on  Displacement  and  Total  Boundary-Layer  Thicknesses 


AEDC>TR.66-37 


Mqj-IO 
ReQj/in.  -396 

^a)^'/^®oo,  i " 

T,^/To  =  0.2  (No  Slip) 

Xqj  =  0.037  in. 

Vg/Ug  Matching  Condition 

All  Calculations  with  TVC 


0.40  0.60  0.80  1.00 


X 


HL 

0 

.0 

0. 198 


b.  On  Stanton  Number 
Fig.  Concluded 


AEDC-TR-66-37 


Experimental  Data 


Sym 

ReQj/ln. 

To,“K 

T^.OK 

□ 

9.38 

2606  ^ 

1033 

-293 

o 

9.37 

1602 

1662 

-293 

A 

10.15 

388 

3020 

-293 

Flagged  Symbols  Uncorrected  Data 

Numerical  Results 

Sym 

^CD 

Re^jj/fn. 

To.  “K 

^w^o 

9.37a 

2869 

1073 

4.95 

0.273 

9.37b 

1632 

1835 

2.76 

0. 160 

10.0 

396 

3034 

2.02 

0.097 

Note:  All  boundary- layer  calculations  include 


Fig.  17  Exporimental  gnd  Theoretically  Predicted  Surface  Pressures 
on  a  Water-Cooled  9-deg  Half-Angle  Cone  at  -  10 


46 


AEDC'TR-«6-37 


To 


u 

+ 

cT 

u 

n 

a 

u 


(2s  +  2s) 


Re^j/in.  ”396 
I "  0*  326 

T^^;ro-o.2(Nosijp) 
>00  ”0.037  In. 
Nitrogen  Injection 


Numerical  Results 
Vg/Ug  Matching  Condition 
□  No  Slip 

A  Slip  and  Temperature  Jump 

Experimental  Results 

AEDC-VKF  Tunnel  L 
- Little  and  Griffith 


^  ±7  percent 


M®-9.37 
Re^/In.  -  1632 


ifi  X  10^,  Ibm/sec 

Fifl,  18  Experimantol  and  Theoretically  Predicted  Zero>Lift  Drag  with  Meet  Transfer 


47 


AEDC-TR-66-37 


MflO  =  10 


ReJ\n.  -22,528 


0.1 


Cd 


0.3 

0.2 

) 

0.1 


T,/ro-a33- 

0.75 

Xqo  -  0. 00066 

Nitrogen  Injection 

Numerical  Results 

Experimental  Results 

□  No  Slip 

AEDC-VKF  Tunnel  "C"  ' 

A  Slip  and  Temp  Jump 

- Little  and  Griffith 

Vg/Ug  Matching  Condition 

IDIls) 

- i  ±7  Percent 

1  1  1  1 

(2)° 

1  1 

0  1.0  2.0  3.0  4.0 

5.0  6.0 

m  X  10^  Ib/sec 

M^-IO 

Reoj/in.  “45,057 

Ma>VReoo,i 

T^o“0-33- 

0.75 

Xgg- 0.00(03  in. 

Nitrogen  Injection 

r  Numerical  Results 

Experimental  Results 

□  No  Slip 

AEDC-VKF  Tunnel  "C" 

A  Slip  and  Temp  Jump 

- Little  and  Griffith 

Vg/Ug  Matching  Condition 

r- 

‘  (lids) 

_ 1 _ 1 _ 1 _ 1 _ 

1  1 

1.0  2.0  3.0  4.0 

m  X  10^  Ib/sec 


5.0  6.0 


Fig.  18  Concigdwd 


48 


AEDC-TR-66-37 


Moo  “10 
Re^j/in.  =396 
-0.326 
T,/To-a2 
Xgj  "0.037  in. 

Nitrogen  Injection 

Numerical  Results  Experimental  Results 

Matching  Condition  Slip  AEDC-VKF Tunnel  L 
d6*/dx  No  - Little  and  Griffith 


£6  d6*/dx  Yes 

o  d6*/dx+ Py^Vy^^r^PgUgfg  No 
u  d6*/dx  +  P^^v^r^PgUgfg  '  Yes 


Fig.  19  Effect*  of  Iteration,  Mas*  Transfer,  Slip,  and  Temperature 
Jump  on  Zero-Lift  Drag  at  •>  10 


49 


AEDC-TR-66-37 


Experimental  Data 


Sym 

Mg 

Re^j/in. 

deg 

To,  "K 

Xqj.  in. 

Tunnel 

Source 

0 

10.15 

388 

9 

3,020 

0.30 

0.630 

0.037 

L 

Little  and  Griffith 

□ 

9.3? 

1.600 

9 

1,662 

0.50 

0.649 

0.0065 

L 

'q 

9.38 

2.600 

9 

1,033 

“0.65 

0.701 

0.0048 

L 

o 

10.00 

23.800 

9 

958 

-0.75 

0.660 

0.00066 

C 

o 

10.00 

42,600 

9 

958 

-0.75 

0.685 

0.00033 

C 

A 

3.93 

1,765 

5 

300 

0.96 

0.821 

0.0032 

Berkeley 

King  and  Talbot 

O 

5.64 

6,200 

5 

300 

0.96 

0.815 

0.0012 

Berkeley 

King  and  Talbot 

Whitfield  and 

Griffith  fl963i 


50 


TABLE  i 

CONDITIONS  USED  IN  THE  EXPERIMENTAL  AND  NUMERICAL  STUDIES 
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P»'* 

Ho  *  lO"’. 

Tar'K.' 

»  lo’, 

-  p«i« 

T_,‘>K 

/>..  *  10*, 

ft/sec 
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3.93 

1765 

0.0334 

0.  3245 

0.  227 
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1.64 

73.4 

1.  04 

2215 

0.821 

5,64 

G200 

0.0679 

0.  3245 

1.775 
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1.64 

40,  7 
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2368 

0.815 

AEDC- 
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L^ 

10. 15 

388 
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3.  930 

18.00 
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0.397 
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0.  124 

8118 
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(10.  00) 
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(3.  395) 

(J6.64) 

(3034) 

0.392 

(0.  122) 

(8041) 
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9.37 

1602 

0.  1056 

2,053 

25. 00 

1662 

0.952 

93.0 
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L 
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(25.62) 

(1835) 
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L, 

9.38 

2606 

0.  0812 

1.200 

19.50 

1033 

0.713 

56.8 

0.565 
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(9.37) 

(2869) 

(19.  70) 
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(0.552) 

(4766) 
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c 
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958 
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45.6 

8.01 

4444 
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(45.  057) 
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(7.57) 

(47.3) 

(7.  44) 
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4.06 
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(994) 

(3.78) 

(3.  72) 

(4526) 

0.675 

^Unless  otherwise  noted  ideal  ^as  data  computed  from  Ha,  p',  and  shown  In  parentheses, 
^deal  gas  data  computed  from  pj,  T^,  and  shown  in  parentheses. 
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TABLE  II 

EXPERIMENTAL  DATA  ON  THE  EFFECTS  OF  MASS  TRANSFER  ON  ZERO-LIFT  DRAG 
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TABLE  III 

EXPERIMENTAL  WALL  PRESSURE  DATA  ON  A  9.DEG  WATER-COOLED  CONE  IN  TUNNEL  L 


Mee 

X 

Pm  ■^Ps« 
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^Pressures  corrected  for  orifice  size  after  Potter,  Klnslow,  and 
Boylan  (Ref.  13). 


TABLE  IV 

EXPERIMENTAL  EQUILIBRIUM  WALL-TO-STAGNATION  TEMPERATURE  DATA  ON 
A  THIN-WALLED  STAINLESS  STEEL  MODEL  IN  TUNNEL  L 
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